Collective excitations in the neutron star inner crust 



O 
(N 

-(— > 
o 

O 

o 






> 

O 

O 



% 



Luc Di Gallo,i'B Micaela Oertel,i'[!| and Michael Urban^'l 

^LUTH, Observatoire de Paris, CNRS, Universite Paris Diderot, 5 place Jules Janssen, 92195 Meudon, France 
^Institut de Physique Nucleaire, CNRS/IN2P3 and Universite Paris Sud 11, 91^06 Orsay, France 

We study the spectrum of collective excitations in the inhomogeneous phases in the neutron star 
inner crust within a superfluid hydrodynamics approach. Our aim is to describe the whole range 
of wavelengths, from the long-wavelength limit which can be described by macroscopic approaches 
and which is crucial for the low-energy part of the spectrum, to wavelengths of the order of the 
dimensions of the Wigner-Seitz cells, corresponding to the modes usually described in microscopic 
calculations. As an application, we will discuss the contribution of these collective modes to the 
specific heat of the "lasagna" phase in comparison with other known contributions. 

PACS numbers: 26.60.Gj 



I. INTRODUCTION 

Neutron stars are fascinating objects, containing mat- 
ter under extreme conditions of temperature, density and 
magnetic field. In order to study these celestial bodies, 
theoretical modeling has to be confronted with observa- 
tions. A prominent observable is the thermal evolution of 
isolated neutron stars. Properties of the crust thereby in- 
fluence the cooling process mainly during the first 50-100 
years, when the crust stays hotter than the core which 
cools down very efficiently via neutrino emission (see e.g. 
[ij). Heat transport in the crust is the key ingredient 
to explain the afterburst relaxation in X-ray transients, 
too [2, 3J. Concerning the models for the thermal relax- 
ation of the crust, the most important microscopic ingre- 
dients are thermal conductivity and heat capacity, and 
to less extent neutrino emissivities. Here, as a first ap- 
plication, we will concentrate on the heat capacity. More 
details about the evaluation of the heat capacity and a 
discussion of the usually considered contributions can be 
found in 4]. In what follows, we will concentrate on 
the particularly interesting case of the neutron star inner 
crust. 

The core of neutron stars is composed most probably 
of homogeneous neutron rich matter, whereas the crust 
contains different inhomogeneous structures. The inner 
crust is thereby characterized by the transition from a lat- 
tice of atomic nuclei in the outer crust to homogeneous 
matter in the core. Ravenhall et al. p] and Hashimoto et 
al. [y] predicted that this transition passes via more and 
more deformed nuclei. Starting from an almost spheri- 
cal shape, they could form rods or slabs immersed in a 
neutron gas at the different densities. These "spaghetti" 
and "lasagna" phases are commonly called the nuclear 
"pasta". At higher densities, even closer to the core, 
other phases such as neutron-gas bubbles inside the dense 
matter ("swiss cheese" phase) etc. are expected. The 
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formation of the different structures strongly depends 
on the relative strength of the nuclear surface energy. 
Coulomb energy and bulk energy, such that it depends 
on the nuclear interaction. This prediction has been con- 
firmed within different models for the nuclear interaction, 
see, e.g., [7HlO|. These evaluations have been performed 
at zero temperature. It is clear that at some critical 
temperature the pasta structures will disappear due to 
thermal excitations. However, this mel ting temperature 
is of the order of several MeV (see e.g. [ll|, [l2|). 

Another point is that in neutron stars older than sev- 
eral minutes, matter becomes superfluid. A first evidence 
for superfluidity in neutron stars has been discussed al- 
ready in 1969 [13, shortly after the discovery of the first 
pulsars, in connection with the observation of "glitches" . 
Since then much effort has been devoted to the ques- 
tion of superfluidity and superconductivity in neutron 
star matter, for the inner crust as well as for the homo- 
geneous core, see for example [8[. There is no consen- 
sus on the exact value of the energy gaps A in the inner 
crust [IJ, ^M, ; but the common agreement is that they are 
of the order of 1 MeV [1^. A pairing gap much larger 
than the temperature strongly suppresses the contribu- 
tion of individual neutrons to the specific heat which is 
thus very much dependent on the pairing strength |4l.ll7|. 
For moderate and strong pairing, the main contributions 
to the heat capacity considered so far in the crust are 
thus electrons and lattice vibrations as well as collective 
excitations of nuclei. However, the superfluid character 
of neutron star matter induces collective excitations of 
the neutron gas, not considered before, which can give 
an important contribution to the heat capacity in cer- 
tain regions, see [18-20]. 

The aim of the present paper is to study these col- 
lective excitations in the inner crust employing a super- 
fluid hydrodynamics approach. Naturally, there is a vast 
literature on hydrodynamics for neutron stars in differ- 
ent contexts, including the effects of superfluidity [211- 
|27| . Most of these models are dedicated to the study of 
macroscopic neutron star properties, whereas our main 
aim is to study the excitation spectrum of the crust on 
much smaller length scales. In spirit this is similar to 
Refs. [1^ [20, [2^ J where hydrodynamic equations are de- 



veloped to study the superfluid Goldstone boson and (lat- 
tice) phonons in the long wavelength limit. However, we 
are interested in shorter wavelengths at which effects of 
the inhomogeneous structure will manifest themselves, 
too. In this sense, our approach is situated in between 
the long wavelength limit and the completely microscopic 
calculations (see e.g. ^3[[2E|33) employing the Wigner- 
Seitz approximation |31| . In the latter case, the wave- 
lengths are limited to the size of the Wigner-Seitz cell, 
because of the imposed boundary conditions which do 
not include the coupling between neighboring cells. 

The paper is organized as follows. In Sec. HIl we de- 
scribe the superfluid hydrodynamics approach. We sum- 
marize the hydrodynamic equations, discuss the bound- 
ary conditions and the microscopic input. In Sec. lIIIi we 
show our first results. For simplicity, we restrict ourselves 
to one-dimensional inhomogeneities (lasagna phase) in 
this exploratory study. We discuss the spectrum of the 
collective modes and their contributions to the specific 
heat. A summary and perspectives of our work are ex- 
posed in Section HVl 

Throughout the article, c, h, and ks denote the speed 
of light, the reduced Planck constant, and the Boltzmann 
constant, respectively. 



II. MODEL 
A. Superfluid hydrodynamics approach 

In this paper, we are interested in temperatures be- 
low ~ 10^ K, which is very small compared to the gap 
energy A. Therefore we can use the zero temperature 
approximation, thus assuming that there are no normal 
fiuids but only superfluids. In this limit, the dynamics of 
a superfluid system with slow temporal and spatial vari- 
ations is completely determined by the dynamics of the 
phase of the order parameter: If the superfluid order pa- 
rameter is written as A(r) = |A(r)|e"^^'22 the superfluid 
velocity is given by Vs — (?i/2to)V</> |32| |. m being the 
nucleon mass. Actually, as pointed out in Ref. [33], the 
phase of the order parameter determines the momentum 
per particle p = niVg and not the fluid velocity v. This 
distinction is important in the context of "entrainment" 
in a system containing protons and neutrons, see below. 

An important length scale is the superfluid coherence 
length, ^0 = Tjwf/ttA [3J|, where vp denotes the Fermi 
velocity. It varies from several fm up to tens of fm for typ- 
ical values of the densities, neutron fractions, and gaps in 
the inner crust. As can be shown by deriving the equa- 
tions of superfluid hydrodynamics from the microscopic 
time-dependent Bogoliubov-dc Gennes (or Hartree-Fock- 
Bogoliubov) equations [35, 36], the hydrodynamic ap- 
proach is valid if the length scale of spatial variations 
is larger than ^Qj and frequencies are small compared 
with A/h. As will be discussed later on, for the concrete 
examples we consider, we are at the limits of validity of 
the approach. However, considering the tremendous dif- 



ficulties to perform completely microscopic calculations 
beyond the Wigner-Seitz approximation, we leave such 
investigations for the future and consider our approach 
sufficient for the moment. 

In addition, wc will neglect the Coulomb interaction of 
the protons. This represents an enormous simplification, 
but at the same time it implies that we cannot correctly 
reproduce the phonons of the Coulomb lattice. In homo- 
geneous matter, too, the Coulomb interaction plays an 
important role for the collective modes, in particular the 
coupling of the proton plasmon mode with the electrons 
as discussed in [33, [S^l , but it is beyond the scope of the 
present paper. Our main focus lies therefore on the dy- 
namics of the neutron gas, which is however coupled to 
the proton dynamics due to the nuclear interaction. 

In principle, the equations of superfluid hydrodynam- 
ics can be derived from the underlying microscopic the- 
ory, as it was done for the case of ultracold trapped 
fermionic atoms in [35,i36|]. Here, we follow the simpler 
way to derive them from local conservation laws. Since 
the fluid velocities and the densities are low enough, we 
will use a non-relativistic formulation. 

The flrst conservation law is neutron and proton num- 
ber conservation [5l|. This results in two continuity 
equations, one for neutrons (a = n) and one for protons 
(a =p), 



dtUa + V • {Ua^a) = 0, 



(1) 



where Ua denotes the particle number density of species 
a. 

The second conservation law, the conservation of mo- 
mentum, results in the Eulcr equations, which can be 
written as 



na (dtPa + ^fj-a) = , 



(2) 



where jla is the rest-frame chemical potential deflned as 
the conjugate momentum with respect to the particle 
density Ua in a variational approach [39|. The explicit 
expression is 



1 2 

A^a = fJ.a + Va • p^ - -UlaV^ , 



(3) 



where /ia is the local chemical potential of species a. Due 
to the interaction between neutrons and protons, ^a de- 
pends on the densities of both species. 

In pure neutron matter, the momentum p„ is simply 
given by p„ = mnVn- However, in a system containing 
neutrons and protons, the two species drag each other 
due to their interaction. In the theory of superfluids, this 
effect is known as entrainment [40]. As a consequence, 
fluid momenta are misaligned with particle velocities. 
The relationship between the velocity and the momen- 
tum can be expressed via the entrainment matrix (also 
called Andreev-Bashkin or mass-density matrix) |33| : 



V- Pb 

^anaVa = > ^ Pab 

/^^ nib 

b^n,p 



(4) 



In practice, at densities which are relevant in the inner 
crust, the non-diagonal elements of p are small [33[, i.e.. 



where 



Pab ~ maTlaSab ■ 



B. Microscopic input 



(5) 



As microscopic input, we need the equation of state, 
i.e., the relation between the densities Ua and the chem- 
ical potentials Paj and the entrainment matrix p. In our 
concrete numerical examples, we will use the results of 
the work by Avancini et al. [lO| for the equilibrium config- 
urations. They evaluate the structure of the pasta phases 
for charge neutral matter in /3 equilibrium using a den- 
sity dependent relativistic mean-field model, the DDHJ 
model (originally called DDEpS) [10, SH, SI , for the nu- 
clear interaction. In order to be consistent, we shall cal- 
culate the chemical potentials pa and the entrainment 
matrix p with the same interaction. For the entrainment 
matrix, we closely follow Gusakov et al. |43;], who gen- 
eralized the determination of the entrainment matrix for 
neutron-proton mixtures based on Landau-Fermi liquid 
theory [4J| to relativistic models. The only modifica- 
tion of the expressions in Ref. |43j we have to perform 
is due to the presence of the isovector-scalar S meson 
in the DDHi5 model, which modifies the Dirac effective 
nucleon mass. In particular, the latter is no longer the 
same for neutrons and protons. Since our hydrodynamic 
equations are formulated non-relativistically, we consider 
only the non-relativistic limit of the entrainment matrix 
{Pab = rnairibC^Yab in the notation of [43]). 



C. Linearization around stationary equilibrium 

In order to proceed we will linearize Eqs. ([Ij and ^ 
around stationary equilibrium. Let us write the different 
quantities as a sum of their equilibrium value and a per- 
turbation, X — ATcq + 5X (in the case of the velocities 
and momenta we will write the perturbation simply as v^ 
and Pq , respectively, since the equilibrium values of these 
quantities are zero). The equations can be simplified a 
lot, since all temporal and spatial derivatives of equi- 
librium quantities vanish (except at phase boundaries, 
which will be treated in the next subsection). Eqs. ([1]) 
and (m then reduce to 



dtSua 



E 

b—n.p 



Pab; 



eq 



marrib 



V-Pfc, 



(6) 



and Eqs. © and ^ become 

dtVa = -^5pa . 



(7) 

We will now express the variation of the densities in 
Eq. ([6]) in terms of the variation of the chemical poten- 
tials, 






(8) 



Jab — 



_ f dUg 

^dpb 



(9) 



Inserting the resulting equation into the divergence of 
Eq. ([7]) one obtains the following system of two coupled 
wave equations for 6pn and 6pp: 



{KJ)ab dfSpb = y^Spa 



b—n.'. 



where K is the inverse of the matrix 



{K-\b^ 



Pab^cq 

marrib 



(10) 



(11) 



The coupling arises from the non-diagonal elements of 
the matrices J and K due to the neutron-proton interac- 
tion. Let us now make the ansatz that the perturbations 
have the form of a plane wave, Spa(j^,t) = Uae~'^'^*'^^^''^ ■ 
Eq. (|10p can then be written as a 2 x 2 eigenvalue problem 



Y, {KJ)abUb = \Ua . 



(12) 



b^n.p 



with u — u)/k denoting the sound velocity. The two 
eigenvalues give two sound velocities which we will label 
u^. Note that the corresponding eigenvectors, U^, do 
not describe pure proton or neutron waves, but combina- 
tions of both. We denote by -I- and — the modes where 
neutrons and protons oscillate in phase and out of phase, 
respectively. 

In the special case of pure neutron matter, there is 
only one mode, which can be obtained from the above 
equations by setting rip = 0. Its sound velocity is given 
by 



(^1^) . (13) 

V m„ Orin ' cq 



D. Boundary conditions 

In our model, we consider the inhomogeneous phases 
in the inner crust as mixed phases where a neutron gas 
(phase 1) coexists with a dense phase (phase 2) contain- 
ing neutrons and protons. However, in order not to have 
to write everything separately for phase 1 and phase 2, 
we will write all equations, unless otherwise stated, for 
the general case that neutrons and protons are present in 
both phases. The equations relevant for phase 1 can eas- 
ily be obtained by considering the special case ripi = 0. 
The fact that both phases coexist implies that in equi- 
librium the chemical potentials and pressures are equal 
in both phases: pa\ — Pa2 and Pi — P2. The descrip- 
tion of the interface between the two phases requires a 
microscopic formalism and is beyond the scope of this 
work. 

In our model, we assume that the hydrodynamic equa- 
tions are valid in both the gas and the dense phase, but 



since they do not say anything about the behavior at 
the interface, they have to be supplemented by appropri- 
ate boundary conditions. The first boundary condition 
arises from the obvious requirement that contact has to 
be maintained at all times at the interface [45|. There- 
fore, the displacement normal to the interface has to be 
continuous and equal for all components (a = n,p) at all 
times. Hence, the velocities normal to the interface must 
satisfy: 



w_L„i(r) = w±pi(r) = w±„2(i") = vj_p2{r) . 



(14) 



The second boundary condition arises from the require- 
ment that the pressure P on both sides of the interface 
must be equal [46| : 

Pi(r)=P2(r). (15) 

If we linearize this condition, it can be written as 

^ nai{r)Sfiai{r) ^ ^ na2{r)Sfia2{r) , (16) 



a^n^p 



a^n,p 



where the index eq after Uai and na2 has been dropped 
for brevity. 

Before applying our model to the neutron-star inner 
crust, let us see whether these boundary conditions give 
reasonable results for collective modes in isolated nuclei. 
For simplicity, we will consider a nucleus with equal num- 
bers of neutrons and protons {N = Z = A/ 2). Within 
the hydrodynamic model, the nucleus is a homogeneous 
sphere with a sharp surface at r = R. The proton and 



neutron densities inside the nucleus are n^ 



no/2. 



where uq — 0.153 fm~^ is the saturation density within 
the DDH5 model. As a first example, we consider the 
isoscalar monopole mode, where neutrons and protons 
oscillate together in radial direction. The solution of 
the wave equation inside the nucleus is 5fi„ = SfXp oc 
jo(ujr/u'^), where ji is a spherical Bessel function and 
M+ = 0.169 c is the sound velocity for the in-phase oscil- 
lation of neutrons and protons. Since protons and neu- 
trons move together, the first boundary condition (ITU) is 
automatically satisfied, while the second one, Eq. ([TB]). 
requires dfi{r = R) = 0. Consequently, the energy of the 
monopole mode is hui = nhu^ /R « 90 MeV/A^/^. 

Another interesting simple case is the isovector giant- 
dipole resonance (GDR) , where neutrons and protons os- 
cillate against each other in z direction. In this case, 
our approach is identical to the Steinwedel- Jensen model 
of the GDR [47]. Again, the solution of the wave 
equation is straight-forward and gives (5/i„ = —SfJ.p oc 
ji{Lur/u~) cos 6, where 9 is the angle between r and 
the z axis, and u~ = 0.233 c is the sound velocity for 
the out-of-phase oscillation of neutrons and protons. In 
this case, the second condition p^ is automatically ful- 
filled, but now the first one, Eq. ([14]), becomes rele- 
vant. Using the Euler equation ([7]), one can show that 
the radial component of the velocity field is proportional 
to frn, = — "Urn OC dS^/dr OC j'i{u!r /u~) COS 9, so that 
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FIG. 1: Diagram representing the slab structure. 



These results for the isoscalar monopole and the isovec- 
tor GDR are quite reasonable, at least for heavy nuclei, 
although their energies are much higher than the pair- 
ing gap A, such that superfluid hydrodynamics should 
strictly speaking not be applicable. The reason is that 
for these particular resonances (contrary to, e.g., the 
quadrupole mode [4^) the Fermi-surface distortion does 
not play any role, so that hydrodynamics works even in 
the normal phase without pairing. We conclude that, at 
least in some cases, the limits of validity of the hydrody- 
namic approach may be interpreted very generously. 



E. Collective modes in a periodic slab structure 

Because of their electric charge, the droplets (or rods, 
or slabs) of the dense phase arrange in a regular peri- 
odic lattice in order to minimize the Coulomb energy. 
Charge neutrality on a macroscopic scale is guaranteed 
by the presence of an almost uniform, strongly degener- 
ate electron gas. The size and form of the structures is 
determined by the interplay of Coulomb energy (favor- 
ing small structures) and surface energy (favoring large 
structures). Since both the Coulomb and the surface en- 
ergy are neglected in our approach, the determination of 
the size and form of the structures in equilibrium is be- 
yond the scope of our work. Instead, we will consider 
the equilibrium geometry as input and calculate the col- 
lective oscillations in this geometry. For simplicity, we 
will restrict ourselves to the simplest geometry which is 
a structure of periodically alternating slabs with differ- 
ent proton and neutron densities as illustrated in Fig. [1] 
(lasagna phase). To be specific, we will consider the slabs 
to be perpendicular to the z axis. 

Our aim is to describe the collective modes of this 
structure. The equilibrium properties of the structure 
itself, i.e., the densities n„i, ripi, n„2 and np2, and the 
slab thicknesses Li and L2, are input parameters which 
we take from Rcf. [lO]. The excitations are then obtained 
by solving in each slab the wave equation (flUj) together 
with the boundary conditions (fTi]) and (fTB)) at the inter- 
faces between neighboring slabs. 

At each phase boundary, the waves will be partially (or 
totally) reflected. It is therefore not sufficient to make a 
plane- wave ansatz in each slab, but one has to consider 
the reflected wave, too. Thus, we can make the following 



ansatz in each slab 






■p^e' 



u: 



(17) 



where U^ denote the normalized eigenvectors of Eq. (IT^ , 
k|| = (fca;,fcy,0) and r|| = (a;,y, 0) are the components of 
k and r parallel to the slab, and the kf have to satisfy 

,2 

(18) 



U±2 



,±2 



'II 



Note that kf can be real or imaginary. The velocities v^ 
can be expressed in terms of the coefficients a* and /3^, 
too. By using the Euler equation ([7]), one finds that for a 
plane wave with wave vector k the velocity field is given 
by 



— V(i^-i),,(5/.,. 

UnUJ ^ ' 

b 



If we define 



V^ = 



b 



{K-%bU^ , 



(19) 



(20) 



the superposition of plane waves according to Eq. (J17p 
gives 






(a^t 



- P^e 



V° 



(21) 
and a similar relation for V||. 

The next step is to determine the coefficients a and 
/3 by matching the solutions in neighboring slabs ac- 
cording to the boundary conditions. If we use indices 
1, 2, 3 in order to indicate the quantities in three con- 
secutive slabs, we have perturbations Sjiai valid for < 
z < Li, 5iia2 for Li < z < L = Li + L2, and 6^ad, for 
L < z < L -\- Li with four unknown amplitudes, a and 
/3^, in each slab. Note that due to the periodicity the 
equilibrium properties of slab 3 are equal to those of slab 
1, but the coefficients a^ and /3^ are in general different 
in slabs 1 and 3. 

Written explicitly, the boundary conditions are 

Vzp2{z = Li) = Vzri2{z = Li) , 
Vznl[z = Li) = Vz„2{z = ii) , 
Vzpl{z = Li) ^ Vzn2{z = Li) , 

^ nalSfJ.al{z = Li) = ^ na2SHa2{z^ Li), (22) 
a—n,p a—n,p 

and four analogous equations relating quantities of slabs 
2 and 3 at z — L. 

It is evident that, in order to satisfy the conditions for 
all r|| , the components k|| must be equal in all three slabs, 
i.e.. 



kill =k||2 = k||3 = q 



^11 



(23) 



So far, the boundary conditions give us eight equations 
for the twelve unknown coefhcients af , • • • , P^ ■ In order 
to close the system of equations, we have to take into 
account translational invariance of the system. This can 
be expressed via the Floquet-Bloch theorem |49j : 



S^iair + -R,t) = e'^■'^S^iair,t) 



(24) 



where q is the Bloch momentum and R = {Rx,Ry,Rz) 
is a vector such that the system is invariant under a shift 
r — !■ r + R. In our case, Ry can be arbitrary, but Rz has 
to be a multiple of the periodicity L, see Fig.[TJ With re- 
spect to R||, the condition (|24)) is automatically satisfied 
due to Eq. ^. But in the case R = (0,0, L), Eq. ^ 
implies in particular 

6naz{x,V,z = L,t) = e"-"^'^ S^iaiix,y,z == 0,t) (25) 

and an analogous relation for the velocity. Inserting this 
into the boundary conditions at the interface between 
slabs 2 and 3 at z = i, we obtain 

Vzp2{z = L) = Vzn2{z = L) , 
e"^'^ Vznl{z = 0) = Vzn2{z = L) , 
e*'^^'^ Vzpl{z = 0) = Vzn2{z = L) , 



gjg.L y^ nal6Hal{z = 0) = ^ na25Ha2{z = L) , 



a—n,p 



a—n,p 



(26) 



i.e., we have now a system of eight equations, Eqs. ([2^ 
and (PSJ) . for eight coefhcients a^^ , . . . , /32 ■ 

This system of equations has a non-trivial solution if 
the determinant of the corresponding 8x8 matrix van- 
ishes. For a given choice of qu and Qz {qz may be limited 
to the first Brillouin zone, i.e., —tt/L < qz < tt/L), this 
gives us an equation for u with an infinite number of 
discrete solutions. 

Note that, as mentioned before, in the case we will 
actually consider, the proton density vanishes in slab 1 
(and 3). In this case, the proton velocity is not defined 
in that slab and we have only two coefficients ai and /3i 
instead of four coefficients ai and /3^ , since in pure neu- 
tron matter there is only one eigenmode instead of two. 
The number of equations is also reduced by two, since 
the third equation of Eqs. (1^ and the third equation 
of Eqs. (Uni) can be removed. We are therefore left with 
a 6 X 6 instead of 8 x 8 problem. In this case, it is in- 
teresting to notice that there are two different types of 
modes: Modes propagating through all slabs, whose en- 
ergies depend on qz , and modes of the dense slabs (slab 2) 
only, whose energies are independent of qz . The latter are 
modes where protons and neutrons oscillate against each 
other such that at the boundaries (z = Li and z — L) 
Vzn, Vzp, and 5P vanish simultaneously (analogous to the 
isovector GDR in an isolated nucleus, discussed at the 
end of the previous subsection) . 

When looking for the solutions for w, one has to be 
careful to retain only physical solutions. It is easy to see 



TABLE I: Properties of the lasagna phase within the model 
by Avancini et al. u}\ studied in our example. The average 
densities of the total system are given by n„ = {L\ / L)nni + 
{L2/L)nn2 etc. Baryon density and proton fraction are de- 
fined as ub — Tin + Up and Yp — Up/rin, respectively. 







slab 1 


slab 2 


total 


L 


(fm) 


9.40 


7.38 


16.78 


n„ 


(fm-«) 


0.0701 


0.0885 


0.0782 


Up 


(fm-«) 





0.0041 


0.0018 


Ub =nn + Up 


(fm-«) 


0.0701 


0.0926 


0.0800 


Yp — Up/UB 







0.0447 


0.0227 


u or u 


(c) 


0.0641 


0.0354 




u~ 


(c) 




0.1369 





that if one of the three wave numbers k. 



k+ 



^z2 



vanishes, i.e., if uj/q\\ equals one of the three sound ve- 
locities ui, U2 , or U2 , the system of equations is solved 
by choosing the corresponding coefficients as a = —(3 
and setting all the other coefficients equal to zero. How- 
ever, this solution implies Sfi = and therefore does not 
correspond to a physical excitation. 



III. RESULTS FOR THE LASAGNA PHASE 
A. Excitation spectrum 

Let us now investigate the resulting excitation spec- 
trum for a specific example. As mentioned before, the 
values for the equilibrium quantities will be taken from 
the work by Avancini et al. [lO|, who have studied the 
structure of pasta phases in a relativistic mean field 
model. Our geometry corresponds to the lasagna phase, 
appearing close to the transition to uniform matter in 
the core, which has been found in ReL [10] in the case of 
zero temperature and /^-equilibrium for baryon number 
densities 0.077 fm^'^ ^ "-s ^ 0.084 fm"-^, in good agree- 
ment with the results by Oyamatsu [7| . For our example 
we have chosen an intermediate density, ub = 0.08 fm~'^. 
The corresponding properties of the two phases 1 and 2 
are listed in Table H] 

With the actual numbers for the densities and the di- 
mensions of the structure, the coherence length for a 
gap of 1 MeV is of the same order of magnitude as 
the size of the layers, i.e. the scale for spatial varia- 
tions. That means that our superfiuid hydrodynamics 
approach touches its limit of validity for this example. 
Strictly speaking, we should also limit ourselves to ener- 
gies which are small compared to A. However, there are 
many cases where the hydrodynamic approach works rea- 
sonably well although its initial assumptions are not ful- 
filled. Examples are the dipole and monopole resonances 
in ordinary nuclei mentioned in the preceding section, 
or the "supergiant resonances" in spherical Wigner-Seitz 
cells used to model the neutron-star inner crust, whose 



excitation energies agree well with an estimate obtained 
from the sound velocity of the hydrodynamic Bogoliubov- 
Anderson mode [30] . 

After this remark of caution, let us discuss the solu- 
tions for the energies w shown in Fig. [2] as functions of 
q = \q\ for three different angles 9 between q and the z 
axis (i.e., Qz = qcosO and gii = qsinO). The left panel 
shows the dispersion relation for waves propagating in 
z-direction, i.e. perpendicular to the interfaces between 
the different slabs. One observes an acoustic branch with 
an approximately linear dispersion law 



Usq 



(27) 



at low energies, and several optical branches with a finite 
energy for q — 0, analogously to phonons branches in a 
crystal. 

Note that within the Wigner-Seitz approximation, 
which is usually employed in microscopic calculations 
pTI . I29I ISOJ I , we would only obtain a discrete spectrum 
corresponding to our spectrum in the case q = 0. The 
reason is that in this approximation the coupling between 
cells is neglected, and thus each cell has the same excita- 
tion spectrum. The degeneracy of the modes in each cell 
is lifted by the coupling between cells, which gives rise 
to a momentum dependent spectrum as obtained in our 
approach. 

The slope of the acoustic branch, i.e., the speed of 
sound, coincides (see dashed line in Fig. [2|) with the usual 
thermodynamic expression for the sound velocity 



1 dP 



m duB 



(28) 



where hb is the average baryon density of the inhomo- 
geneous phase. To evaluate this derivative, we squeeze 
or expand our unit cell of length L by a small amount 
6L — 6L1 + 6L2. From the requirement 6 Pi = 6P2 = SP 
we can determine SLi and SL2 and thus SP. The final 
result can be written in a compact form as 
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where we have defined for each phase i = 1 , 2 



1 dP, 



m duBi 



Y^, 



(29) 



(30) 



Note that Ugi is identical to the sound velocity ui [cL 
Eq. (jlSpj. whereas Us2 is different from the two sound 



velocities Uj . 

This linear branch corresponds rou ghly to the long 
wavelength limit discussed in Ref. [10, [2d| although, of 
course, the numerical value of the sound speed is not the 
same because we neglect elastic effects of the proton lat- 
tice due to Coulomb interaction. At higher wave vectors 
q, there are deviations from the linear behavior related 
to the inhomogeneous structure. At these energies the 
long wavelength limit is no longer valid. 
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FIG. 2: Dispersion relation of the modes propagating along the z-axis {9 = 0, left), at an angle of 45° (6 = 7r/4, center), and in 
the plane parallel to the slabs {6 = n/2, right). The dashed line in the left panel corresponds to the approximation Eq. (|27p . 



The central and right panels of Fig. [2] show the excita- 
tion spectrum for different values of the angle, 6 — 7r/4 
and 7r/2, respectively. One observes that the slope of the 
acoustic branch discussed before changes: in the present 
example, Us increases from 0.072 c in the case 6* = to 
0.085 c in the case 9 = 7r/2. The reason is that the wave, 
which is perfectly longitudinal (v || q) in the case 9 = 0, 
becomes more complicated in the case 6* 7^ and the 
nucleons oscillate now in both longitudinal and trans- 
verse directions. But the most important consequence of 
non-zero angle 9 is the appearance of a second acoustic 
branch, whose slope is strongly angle dependent. In fact, 
if one writes the energy of this new branch as 



volume, is defined as 



l^ll = u'^qsm9, 



(31) 



the "two-dimensional sound velocity" u'^ defined by this 
equation depends only weakly on q^ and qy : in the 
present example, u'^ varies between 0.04 c for q^ ^ ^n 
and 0.046 c for gn -C g^. A detailed analysis of the solu- 
tions for the coefficients a and /3 corresponding to this 
branch shows that in this mode, the protons and neu- 
trons oscillate practically only in the direction parallel to 
the slabs (i.e., Vz ~ 0), and the motion takes essentially 
place in the dense phase. 



B. Application to specific heat 

We are interested here in the contribution of the above 
discussed excitation modes to the specific heat. The spe- 
cific heat, the heat capacity for constant volume per unit 



CAT) = § 



(32) 



where e denotes the energy density. The contribution of 
the collective modes can be calculated as follows: 



in-CtJM^M^ 



1 



g?i"i(q)/fc£ 



1 



(33) 
Note that we suppose here that the energies Wi(q) de- 
pend only very weakly on temperature such that it is 
justified to neglect their temperature dependence. This 
should be a good approximation as long as the tempera- 
ture stays well below the value of the energy gap and we 
therefore have no significant contribution from a normal 
fluid. Another type of temperature dependence could 
arise from a change in the structure of the pasta phases. 
At the temperatures considered here, however, we do 
not expect a significant effect either since the structure 
starts to be modified considerably only at higher temper- 
atures [ni[ii. 

In Fig. [3] we show the different contributions to the 
specific heat in the density range where the model by 
Avancini et al. [101 predicts the lasagna phase, for a typ- 
ical temperature of 10^ K. Besides the contribution of 
the collective modes (solid line), we display for compari- 
son the contribution of the electrons (dashed line) , which 
are considered as a practically uniform ultra-relativistic 
{pe 3> meC^) ideal Fermi gas with number density rie = 
Up. At low temperature, the electron gas is strongly de- 
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FIG. 3: Different contributions to the specific lieat in tlie 
density range where one expects to find the lasagna phase, 
for T = 10® K. Solid: collective modes, dashed: electrons, 
dotted: neutron quasiparticles (from Ref. [13] )• Concerning 
the conversion between astrophysical and nuclear units, note 
that 10® K = 86.17A;-^ keV and 10^** erg K'^ cm"^ = 7.246 x 
lO'^fes fm-^ 



generate and its contribution to the specific heat reads 



[hcf 



(34) 



The importance of the collective modes becomes clear 
if one considers the contribution of the gapped neutron 
quasiparticles (dotted curve), taken from Ref. 1131: In 
the absence of collective modes, an excitation of the neu- 
tron gas requires the breaking of Cooper pairs, which is 
suppressed by a factor of the order of ^-^I^bT ^ Even 
in the case of weak pairing, at the present temperature, 
this contribution is suppressed by approximately one or- 
der of magnitude with respect to the contribution of the 
collective modes. 

In Fig. Ul we show the temperature dependence of the 
specific heat corresponding to the intermediate-density 
case discussed in Sec. IIII Al (solid line). For comparison, 
we again display the specific heat due to the electrons 
(dashed line). Due to its linear temperature dependence, 
Eq. (j34[) . the electron contribution is always dominant 
at low temperature, but at higher temperature, the con- 
tribution of the collective modes is comparable or even 
larger than the electron contribution. 

At the low temperatures considered here, which are 
well below the energy of the first optical branch, the con- 
tribution of the collective modes to the specific heat is 
completely dominated by the two linear branches dis- 
cussed in the preceding subsection. As is well known 
|50j . the specific heat due to an acoustic branch with a 
linear dispersion relation, Eq. (j27|) . reads 



^ EE hT^ . 



(35) 



In the present case, however, we have seen that there is 
in addition a "two-dimensional" branch which propagates 



FIG. 4: Temperature dependence of the contribution of the 
collective modes (solid line) and of the electrons (dashes) to 
the specific heat for the the example studied in the preceding 



subsection (see Table |l]). The approximate formula , aT"^ 
bT'^, see Eqs. p5|) and (|36|l . is shown as a dotted line. 
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only parallel to the slabs and whose dispersion relation 
is approximately given by Eq. (j3ip . The contribution of 
such a mode to the specific heat is readily shown to be 



3C(3)fc|r^ 



aT^ 



(36) 



where ( is the Riemann zeta function [C(3) = 1.202 . . .]. 
Due to its quadratic temperature dependence, this is the 
next dominant contribution at low temperatures after the 
electrons. The result of the simple formula aT^ + bT^, 
where a and b have been calculated with the average 
values Us = 0.078 c and u'^ — 0.042 c, is shown as the 
dotted line in Fig. 01 Up to the temperatures considered 
here, it fits reasonably well the full calculation. 



IV. SUMMARY 

In this paper, we have presented a formalism of super- 
fiuid hydrodynamics to treat density- wave propagation in 
inhomogeneous pasta-like nuclear structures which ap- 
pear in the inner crust of neutron stars. To account 
for the periodicity of the structure, we incorporate the 
Floquet-Bloch boundary conditions. The idea is some- 
where in between the approaches of Refs. [1^, [20j, con- 
sidering only the long- wavelength limit, averaging over 
the microscopic details of the structure, and microscopic 
calculations of the crust within the Wigner-Seitz approx- 
imation [2^, [so] , valid for wavelengths smaller than the 
radius of the Wigner-Seitz cell. Concerning the micro- 
scopic input for the nuclear equation of state and the ge- 
ometry of the structure, we followed the work by Avancini 
et al. 10]. 

Within this approach, we have calculated the excita- 
tion spectrum of a periodic structure of parallel slabs, 
the lasagna phase. We have shown that the structure 



can indeed induce non-negligible effects on the excita- 
tion spectrum. In particular, we found that the sound 
velocity of the usual acoustic mode depends on the di- 
rection of the propagation, and, more surprisingly, that 
there is a second acoustic mode whose dispersion relation 
is almost independent of qz- In addition, we found dif- 
ferent optical branches, similar to the phonon spectrum 
of ordinary crystals. 

We have calculated the specific heat corresponding to 
this excitation spectrum and found that its contribution 
is much more important than that of individual neutrons, 
which is strongly suppressed due to the superfluid gap 
A. At temperatures relevant for neutron stars, the main 
contributions to the specific heat come from the electrons 
and from the acoustic collective modes. The latter cannot 
be obtained within the Wigner-Seitz approximation. Due 
to the curious sound mode whose energy is independent 
of (7z, the specific heat due to the collective modes goes 
like T^ instead of T^. [With the same arguments, one 
would predict that in a rod structure (spaghetti phase), 
the specific heat should be linear in T.] However, it is 
not clear whether this feature survives when the Coulomb 
interaction, which has been neglected here, will be taken 
into account. 

Of course, in order to treat the complex geom- 
etry, we were obliged to make a couple of ap- 
proximations. Contrary to the microscopic ap- 
proaches based on the ^Quasiparticle- Random-Phase- 
Approximation (QRPA) [23, [sol, we rely on the assump- 
tion that the modes can be described hydrodynamically, 
which implies in particular that the local neutron and 
proton Fermi surfaces stay spherical at all times. This 
assumption is justified if all spatial variations are slow 
compared to the superfiuid coherence length and the tem- 
poral variations are slow compared to the superfiuid gap. 



Both assumptions are not very well fulfilled. However, 
we have cited examples where hydrodynamics gives rea- 
sonable answers even beyond these very restrictive lim- 
its, and we believe that the results are at least qualita- 
tively correct. The most serious limitation of the present 
work is probably that the Coulomb interaction has been 
neglected. The Coulomb interaction between the pro- 
tons results in an additional coupling between neighbor- 
ing cells, which can have important consequences for the 
excitation spectrum. In the approaches of Refs. |19ll20|. 
it was accounted for by including the elasticity of the 
Coulomb lattice. In our more microscopic approach, the 
Coulomb potential would have to be included from the 
beginning into the proton chemical potential iip{r,t) in 
the Euler equation ([2]) . This is a difficult task which will 
be left for future studies. 

It has to be stressed that the contribution of the col- 
lective modes studied here is potentially more important 
than other contributions, notably the contribution from 
individual neutrons. Therefore it is interesting to pursue 
their investigation and to include the additional contribu- 
tion to the specific heat in studies of neutron star thermal 
evolution. 
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